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ABSTRACT 

A theory of the lo-independent decameter radiation is developed. 
The radiation results from excitation of the electromagnetic loss-cone 
instability by keV electrons, stably trapped near L = 6. The radiation 
is excited in Band III of the extraordinary mode. When the effects of 
refraction are estimated, it is shown that above 10 MHz radiation is 
beamed into the equatorial plane in a wide, but thin, conical sheet 
{'if S£ 80°). When the instability analysis is coupled with one of the 
octupole models of the Jovian magnetic field, the maximum convec- 
tive growth of the instability occurs in the directions of the non-lo A, 
B and C sources. The shape of the peak radio flux frequency spectrum 
is found to be a consequence of the loss cone shape of the electron 
distribution function. 
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AN EMISSION MECHANISM FOR THE lO-lNDEPENDENT 
JOVIAN DECAMETER RADIATION 


I. INTRODUCTION 

In the preceding paper (Goldstein, Eviatar, and Thieman 1 978 - hereafter Paper II) 
we introdu^d a model for the lo-independent Jovian decametric emission (DAM) predica- 
ted on four basic assumptions: the radiation is produced by stably trapped electrons; it 
is emitted in the extraordinary mode near the election gyrofrequency, l£2el; it is beamed 
into a thin conical sheet; and finally, the occurrence probability of DAM noise storms is 
associated with the magnitude of the tilt of tl^ emission cone toward the equatorial plane. 
This last assumption accounted well for the northern hemisphere locations of the Early 
(B) and Main (A) sources but did not satisfactorily account for the location of the source 
in the southern hemisphere - the Third source (C). That paper was based solely on geo- 
metrical constructions coupled with use of the high order, non-dipolar models of the 
magnetic field developed following Pioneer 1 1 encounter with Jupiter (Aciina and Ness 
1976, Smith, Davis, and Jones 1976). 

The results of Paper II demonstrated that a detailed knowledge of the Jovian magne- 
tic field in the source regions does provide a key to understanding details of the geometric 
beaming pattern of the lo-independent DAM. However, the model itself was only a 
framework within which a more complete theory of the decameter radiation could be de- 
veloped. Ideally, such a theory would provide a Physical basis for some of the assump- 
tions used in Paper II; that is the purpose of this communication. This work is motivated 
by the results of Paper II; and fits within the general framework of that model. The 
theory is developed from the single assumption that the inner Jovian magnetosphere (L 
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between 2 and 8) contains a significant number of non-relativistic trapped electrons. This 
assumption is indirectly supported by electron measurements made on Pioneers 10 and 
11. For although instrumentation was not available to detect keV electrons, electrons 
with eneigies above several hundred keV were detected throughout the Jovian magneto- 
sphere by Fillius (1976), Van Allen (1976) and McDonald and Trainor (1976). Further- 
more, Van Allen (1976) has reported several electron spectra taken between 3 and 20 R2^ 
which extend down to 60 keV. Thus we feel justified in assuming the existence of trap- 
ped electrons with keV energies in the inner magnetosphere. Our analysis (and notation) 
is closely tied to that of Paper II in that extensive use is again made of the octupole 
models (O4 and JPL) of Jupiter’s magnetic field. 

The excitation mechanism upon which this paper is based is the electromagnetic 
loss-cone instability. We first mentioned this instability in connection with the lo- 
independent DAM in an earlier paper (Goldstein and Eviatar 1972 — hereafter Paper I). 

At that time we assumed that the equatorial surface magnetic field could reach 20G, 
which led to the conclusion that the lo-independent radiation a>uld be generated by elec- 
trons in the magnetic equatorial plane - where the energetic electron densities were 
greater than elsewhere along a magnetic field line. At the small L values considered 
(L 2), the pitch angle anisotropy was large enough to produce substantial growth rates 
of cyclotron waves. 

We now know that the maximum surface magnetic field on Jupiter is about 14G 
(Smith et al„ 1976; Acuna and Ness, 1976), reached only in the northern hemisphere. 

In the south, the fields do not exceed 10.5G, corresponding to an electron gyrofrequency 
of only 29 MHz. In view of these observations it is clear that if the non-Io DAM is 



ampiined cyclotron radiation, it cannot be produced near the equatorial plane, but must 
be ^nerated at fairly high latitudes and low altitudes (Paper II). 


The Pioneer 10 and 1 1 observations have thus prompted us to reexamine the analy- 
sis presented in Paper I, and to propose a somewhat more detailed description of the 
emission mechanism. This paper also provides an opportunity to clarify some confusion 
that has arisen in the literature as to the propagation mode of the unstable electromag- 
netic wave which we had investipted, i.e. whether or not it must tunnel through a stop- 
zone as it propagates away from the source region. 

We pointed out in Paper II that the highest magnetic fields reached by stably trapped 
electrons correspo’.id to electron gyrofrequencies of no more than 29 MHz — close to the 
maximum observed frequency of the non-lo DAM (Wilson, Warwick, and Libby 1968). 
This is consistent with excitation of DAM in the R-X mode very close to either the up- 
per hybrid frequency, c*>uh cutoff frequency, ccr , where (Stix 1962) 


CU?,u “ 

^UH e '"pe 


and 


U.R = '/2in.i(i + (1 + (1) 

where = V4irNevM is the electron plasma frequency. Because the 
everywhere in the inner magnetosphere and ionosphere, lOgl, ojr and nearly 

equal. 

The analysis of the electromagnetic loss-cone instability in Paper I indicated that 
rapid exponential growth of the linearly unstable waves occurred at frequencies, w, just 
below wjjH and just above (Paper I, Figure 2). Waves excited at frequencies below 
couK propagate in what is called the “Z-mode”, or as we prefer, Band II, while the free 
propaption zone above cor will be denoted Band III. 
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It was pointed out in Paper 1 that the growth rates of the instability were compara* 
ble for waves excited in both Bands II and III. However, waves in Band II cannot escape 
the Jovian magnetosphere without encountering the stop zone that is reached when the 
local upper hybrid frequency equals the frfK)uency of the outward propagating wave. For 
situations in which wm] st , a wave in Band II may be able to ‘*tunner* throu^ the 
stop zone (Budden 1955), a possibility that was discussed in Paper I and more recently 
by Smith (1976). Waves excited in Band 111 have frequences above the cutoff, and be- 
cause both the density and field decrease with increaang radius, these waves can freely 
escape the Jovian environment. The loss-cone instability was shown in Paper I to excite 
waves in both Band*; II and III. We emphasize this point because in several sub^uent 
papers this fact was overlooked (Melrose 1976, Smith 1976) and the electromagnetic loss- 
cone instability was described as exciting only Band II emission. In this paper, to avoid 
any further confusion, we will not discuss excitation of Band II emission. The growth 
rates for emission into the (free-escape) Band III are adequate to explain the lo-indepen- 
dent DAM and tunneling appears to be unnecessary, though it may well make a signifi- 
cant additional contribution. 

In the following section we first derive the growth rate for the loss-cone instability 
at any point along a field line. The trapped electrons are assumed to have a distribution 
function that is composed of a pov er-law in velocity together with a loss-cone in pitch 
angle. Smith’s (1976) contention that waves cannot be excited in Band 111 through reso- 
nant interaction is shown to b<* erroneous for very underdense plasmas. The electrons 
and waves arc able to satisfy the appropriate resonance condition in spite of the fact that 
the wave phase velocity exceeds the speed of light. The growth rates are calcula.ed a* 
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several frequencies throughout the Jovian magnetosphere, utilizing both the O 4 and JPL 
spherical harmonic (octupole) models for the magnetic field. 

The instability exdtes electromagnetic radiation above IQ, | in the right circularly 
polarized extraordinary mode (R**X). The effects of refraction on this mode can be esti* 
mated in a straightforward way, as can the total amplification along a ray path. In the 
absence of nonlinear effects, wav» will be amplified until the local value of no 
lon^r permits the wave>particle resonance (xmdition to be satisfied. At frequencies 
above 10 MHz, we show that the radiation pattern is confined to a thin conical sheet 
with half-angle, as 80*85®. At 10 MHz and below, 'P can be as small as 45-50®. 

At frequencies below 20 MHz, the model of Paper II suggested that radiation beam- 
ed into a wide conical sheet could intersect the equatorial plane at all central meridian 
longitudes. Nonetheless, the three decameter sources at those frequencies are still con- 
fined to relatively narrow bands of central meridian longitude (CML). The location of 
the two northern hemisphere sources, non-Io A and B were associated with places where 
the angle, t, between the axis of the emission cone and the equatorial plane, was a mini- 
mum. A »milar suggestion apropos of the lo cx>ntrol!ed sources had been made by Dulk 
(1965) and Schatten and Ness (1971). This idea, which apparently worked so well for 
sources A and B, gave rather puzzling results when applied in the southern hemisphere. 
Although the minimum value of t did occur when one side of the emission cone pointed 
toward the observed CML of non-Io C, it was the side of the cone pointing toward lower 
longitudes. In contrast, lo-C is observed when the phase of lo is near 240®, from which 
one infers that the radiation is beamed into the high longitude side of the emission cone. 
Thus, the model of I^per II, while correctly locating the northern hemisphere sources, 
was not consistent with the inferred location of io-C. 
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We show in § III that the loss-cone instability growth rates in the northern hemi- 
sph^e are kige^ exactly where i is a mmimum, thus providing a natural explanation for 
the fact that DAM is (X>nrin^ to only those central meridian longitudes where t is 
small - elsewhere the growth rates are n^Ugible. By contrast, in the southern hern- 
sphere lai^e growth rates occur near the location of large values of i - the resuming emis' 
sion pattern then beams radiation into the Source C longitude range, but now from the 
same ade of the emi^ion cone as had been previously deduced for lo-C. These results 
obtain only when the O4 magnetic field model is used. The reasons for this are discussed 
below. 

Several additional features of the observations are discussed in § III, including the 
shape of the radio intensity versus frequency relation, and the size of the source region. 
We also present evidence suggesting that the non-Io DAM is produced along the lo flux 
tube (L = 6). 

II. THEORY 

We begin this section with a brief review of the theory of the electromagnetic loss- 
cone instability. The reader is referred to Paper I for more detail. We will be interested, 
as before, in propagation nearly orthogonal to the background magnetic field, and so we 
consider wave vectors k such that k^/kn = tan $ 1 , where |! and i are defined with 

respect to B. When relativistic effects and proton dynamics are ignored, when the bulk 
of the plasma has a low but finite temperature, and when the number density of the 
thermal component is much greater than that of the high-energy component, the growth 
rate of electromagnetic waves is given by (Paper I) 
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are ordinary Bessel functions of integral order, and f^ is the velocity spac^ distribution 
function of the energetic electrons normalized so that 

/d’pf„(p)-f (3) 

where is the latio of the number density of the high*^nergy electrons to that of the 
thermal electrons. The rest of the notatk>n follows that of Papers 1 and II. 

The frequency oi is related to k through the Appleton^Hartree dispenion relation 
(Stix 1962, p. 380 

2vc>)^(l -p) 

n2 = 1 (4) 

2w^(l -i^)- nj sin^ f - lOglA 
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where 

A - {nj sin^ i + 4u>^U cos^ (4) 

and 

P a 

One should note tiiat v e wiU be interested in computing the group velocity, Vg, of the 
unstable waves. This involves taking derivatives with respect to and so we cannot make 
the usual quasi-trans^rse approximation in which A a ^^1 

Before we discuss the f*:nctional form of ) and the evaluation of the growth 
rate, we want to address tht point raised by Smith (1976) that a cyclotron resonance 
instability cannot excite waves in Band 111 (w > In Band III, ce/k > c (cf. Eq. 4); 

’ on the other hand, to excite the instability V|jO = (to + ne)/k|j must be less than c. 
Even with k | <SC k^ the resonance condition c^ be sa*isfied if w is sufficiently close to 
In^l. In general this is difficult because to must also exceed to^. However, where tOp^ <5 
Ifl^l the resonance velocity is approximately given by 

V||® ^ ^ ^ j tan { (5) 

The source region of DAM will lie above the Jovian ionosphere where the electron densit; 
is completely u»->known. However, N will certainly be much less than the maximum iono^ 
spheric density (-'lO^cm'*^) deduced from the Pioneer 1 1 radio occultation measurement 
(Fjeldbo et al. 1976). By 3000 km altitude above the cloud tops the analysis of Fjeldbo 
et al. indicates that N < lO^cm”*^. We consider below the dependence of the loss-conc 
instability growth rates on values of N ranging from 10 to lO^cm"^. As long as N is « 
iC^cm*"^, excitation of ti. ■ extraordinary mode in Band III is not a problem, although 
the bandwidth in a> over which y is appreciable is very narrow. Note that tliroughout 
most of the earth’s magnetosphere ccp,. is the same oruer as 1X2^1, precluding excitation of 
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cyclotron waves directly into Band III. However, as noted by Smith (1976), in the hi^- 
htitude re^n above the plasmasph^ <K 1. This is the source location of the 

Terre^rkl Kitometric Radiation (TKR) (Gumett 1974, Kaiser and Stone 1975), 

As in l^per I we choose the electron distribution to be a loss-cone but with a power 
law velodty dependence. The results of our analysis are not sensitive to the particular 
functional form chosen for the loss cone and, as before, we assume that relativistic effects 
can be neglect^. Thus 

fo = C(1 + p^/pj)’“ fin (sin o/sin Oq) (Oq < a < ^ -Oo) 

-0 (0<a<oco,a-OQ<a<a’) (6) 

where a is the pitdi angle with respect to the magnetic Held, Og the loss-cone angle, pj^ 
a parameter that determines the momentum at which this distribution merges into the 
thermal background, and C is given by 

2"’-* (m-D! f 

C = 

ff^pjlfinctn (-^) - cosotgKZm-S)!! 

In Equation (6), m is a positive integer greater than one. 

We pointed out above, based on the Pioneer 1 1 models of the Jovian magnetic field, 
DAM can be produced at frequencies close to only if the electrons reach the very 
low altitudes and fairly high latitudes at which the gyrofrequency is in the range of the 
non-Io emission frequencies (1 - 28 MHz). Using a model of the magnetic field, the 
mirror altitude and mirror field intensity, the equatorial loss-cone angle, and ff, at the 
magnetic equator can be determined. Liouville’s theorem then prescribes f^ at any mag- 
netic latitude on that field line. 

The problem of mapping trapped particle distributions along field lines has been 
described in detail by Roederer (1970). If we use a subsc*‘ipt “E” to denote quantities 
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evaluated at the magnetic equator, f^(a) at some point along the field line is given by 



(aL< a < 


= 0 (0 < 0( < of^, w-Ol < a < ff) (7) 

where is defined by 

sm “ ( — ) sin 0£g 
Be 

In Equation ( 7 ) B£ is the magnetic field at the magn tic equatorial plane, while B denotes 
the intensity of the field at some point along the field line. At that point, a is the pitch- 
angle and is the local loss-cone angle in the sense that there are no particles at that 
point with a < . 

Once a specific representation is chosen for the magnetic field, the growth rate can 
I valiiated from Equation ( 2 ), using Equations ( 3 ) and ( 4 ) together with Equations (6) 
and ( 7 ), 

The Jovian magnetosphere does not provide a homogeneous environment in which a 
wave can grow at the rate 7. Because of the strong magnetic field gradients, waves ini- 
tially amplified at the rate 7 will soon propagate into regions where the magnetic field 
and hence 7, have changed. Thus 7 will be a function of Although the field gradients 
are large on scales of a Jovian radius, they are generally negligible on a scale of the inverse 
of the maximum unstable wavenumber, so that the variation in 7 with changing can be 
treated fy taking the instability as convective, with an unstable wavenumber K = 7/vg, 
where = dcc/Dk is the group velocity of the wave. To evaluate K,Vg must be known 
as a function of oj and 
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The magnitude of Vg is given by (Stix 1962» p 520 
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where 


-1 

tan 5 = — 9n/d{ 
n 

The angle 5 is the angle between Vg and w/k (Figure 1). From the dispersion relation 
(Eq. 4) and the definition of Vg, one finds 
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(2(1-v)- 5?| sin2{ + 2sf |2 
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< a 3‘ 


and 


tan 5 = 


2vw^(l-i')ng sin^(X)s^ fig sin2|-2cu*(l-p)^ 

[I ] 


n“ 1 2cu( 1 sin“| + 12, A] ^ 


(9) 


n,A 


For parameters thought typical of the Jovian magnetosphere, 6 is generally small 1-2*^) 
and negative, but zero at ^ = 0 or ir/2. 

As a wave begins to propagate and grow, it is likely to be refracted because of the 
strong gradients in B. While it is not our intention to perform ray-tracing calculations, 
we have found that it is quite possible to estimate both where refraction of the extraordi- 
nary mode is important, and by how much a wave is likely to be refracted. Such an 
estimate is crucial for an understanding of just why the decameter radiation is beamed 
into a wide and thin conical sheet with half angle ^ 80° and A'l' s 5°. 
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The theory of \^^ve refraction in a cold plasma has been treated by Stix (1962, 
p. 57-59) and more recently by Smith (1973). We will use only those elements of their 
analyses most germane for our requirements. 

Let D (k, w,^) - 0 be a slow varying function of£. D can be either the dispersion 
relation or some function of the dispersion relation (Smith uses D » co^n^ - c^k‘). 

Then 

dk 3D/3r 

^ (10) 
dt dD/3cj 

For the extraordinary mode 

dD an" 

— = 2co(n-+X ) (11) 

3w ax 


dD 3n" dn" d 

— af 1-2(X + e —) - Sn B1 (12) 

ar ax ae ar 


where 

X =: 

€ = oilJnl 

The eti nation for dD/ajr is valid in regimes where e « 1 and gradients in B are larger 
than gradients in N (Smith 1973). From Equations (10) - (12) we have 


with 
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In the Jovian magnetosphere, the dominant component of the gradient B is atong g. Thus, 

Akii H AB 

At a < (14) 

k n(Vg/c) B 

When u^g the octupole models of the magnetic field n, Vg, and AB/B can be 
calculated as a wave is convectively amplified. Equation (14) permits estimation of the 
amount of refraction likely during the wave’s growth. 

In the following section we utilize the octupole models of the field together with 
the instability analysis to apply the theory to Jupiter. 

III. APPLICATION TO JUPITER 

In Paper II, two of the octupole models of the Jovian magnetic field were discussed 
in some detail. One was the 23 coefficient model constructed by Smith et al. (1976) 
(denoted JPL) which included three internal (up to octupole) and one external (dipole) 
term. The second one was the “O 4 ” model developed by Acuna and Ness (1976) which 
also had three internal terms, but no external ones. These models differ primarily in the 
placement of the high field region in the southern hemisphere. This difference necessi- 
tates using both field models in the instability analysis. 

Before immersing ourselves in detailed calculations of y and K at various latitudes, 
longitudes and radii, it is informative to compute y at one position using various choices for 
the several parameters (^, N, p^^ = Mu^^, m and f). 

As an example, consider the growth rate near 20 MHz at System HI Longitude Xjjj = 
198®, © = 55® north latitude and r = 1.02 R 2 ^ (where R 21 is one Jovian radius). This is 
a point on the L = 6 flux-tube. A plot of y/IOgl versus cj is shown in Figure 2 for Band 
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Ill using { = 80“, N = 100 cm“^, * 2 x 10'^ c (corresponding to 100 eV), m = 3 and 

^ - 0.1. As mentioned above, y is relatively large over only a narrow range of frequencies 
just above . 

We define the maximum value of 7 (a)) as . K (o>) is similarly sharply peaked in 
Ci}, and we define as its maximum value. (The frequency which defines Kj^ is nearly 
equal to that which defines 7 j^ .) In Figure 3, Kj^ is plotted against 5 at the same loca- 
tion. Separate curves are shown for m = 2, 3 and 4, As in Figure 2, u^ = 2 x 10"^ c, 

N = 100 cm"^, and the frequency is 20 MHz. 

In Figure 3, one can begin to understand both why the conical sheet into which the 
radiation is beamed is thin and why 'i' is large. At ^ = 80®, Kj^ has dropped to little 
more than 1/10 of its maximum value near | = 89®. This strong dependence of Km on J 
is primarily due to the equally strong dependence of Vg on which completely dominates 
the tendency of 7 j^ to peak at smaller values of $ (Figure 4). We have also investigated 
the variation of Kj^ with Uj^, N and At 20 MHz, with f = 89®, Kj^ decreases by a 
factor of 3 as u^^ changes from 3 x 10® cm/s to 1.5 x 10^ cm/s. The variation with N is 
somewhat more complicated. In general Kj^ increases with increasing N, from 1.4 x 10"^ 
cm“^ with N = 10 cm"^, to 1.1 x 10"^ cm"^ with N = 100 cm“®. As N increases to 
10 ^ cm‘^, decreases to 3.3 x lO'^ cm“* because is increasing, making the reso- 

nance condition more difficult to satisfy. In fact, for densities in the source region much 
in excess of 10^ cm“^, the loss-cone instability cannot be excited. From Equation ( 6 ) 
one finds that 7 a f. For the remainder of the discussion we will take N = 100 cm“®, 

Uj^ = 2 X 10 '^ c, f = 0.1 and m = 3 . 
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(a) Refraction 

Before discussing the variation of Km with longitude, it is necessary to estimate the 
effects of refraction on the convectively ^prowing wave. As before, we adopt the notation 
and definitions used in I^per 11. Note that longitudes on Jupiter, or in the Jovian mag- 
netosphere, are denoted Xnp where as always, the subscript refers to System III (1965). 
The observed decameter source locations always refer to the central meridian longitude 
(CML) of Jupiter at the time of observation. As in Paper 11, a distinction is made be- 
tween the computed direction in the equatorial plane toward which radiation is beamed 
in the conical sheet and the longitude of the physical location of the emitting region 
(Xiii). The former is denoted Aj|j(±#i) (± designating the side of the cone into which 
the radiation is beamed). DAM beamed into the direction Am will be observable at Earth 
when the Jovian CML equals Ajij. In the discussion below ^ refers to the emission cone 
angle as it might be seen by a remote observer - after the effects of refraction have taken 
their toll. The angle ^ refers to the cone angle in the source region , and will be us«J, 
when needed, to single out the initial value of ( at the apex of the emission cone before 
any refraction has occurred. Thus, a ray might start off with = 89°, be refracted to 
^ = 85° in the middle of the source region, and finally emerge with ^ = 80°. 

In Section II we derived a formula (Eq. 14) useful for estimating the amount of re- 
fraction experienced by a wave propagating through source region. Equation (14) then 
proNides an indication of the ray path. Along the ray path, beginning at the apex of the 
emission cone, the quantity 

s 

A^(S) = 2 fdsK^a) (15) 
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can be computed, at least approximately. In Equation (15), + and - signs refer to the 
two dire^ions toward which the emission cone intersects the ^uatorial plane (i.e. A|||(+m) 
<wr Ajij(^)). At S, B(s) has changed from its value at s = 0 (the apex of the cone) so 
that t(s = S) a 0. A* is roughly proportional to the logarithm of the intensity. In esti- 
mating A^ we have not included any nonlinear saturation processes. In fact, for some 
choices of {diameters, it is not clear that nonlinear processes are r^uired because the 
gradients of the magnetic Held »em adequate to limit the convective growth to that de- 
manded by the observations. We will return to thi*^ point below. 

Because 5 is generally a snmll angle, the important parameter in estimating the refrac- 
tion is L. We have seen in Figure 3 that for « 80'*, is small. In addition v. = c 
and not only is A* small, but refraction is negligible. U«ng Equations (14) and (15) we 
found that when ranged from 85® to 89®, refraction would result in the ray leaving the 
source region with ^ between 80-85®. Re<^ll from Paper II that s 80® was required 
to properly locate the non-Io A and B sources in the northern hemisphere. 

So far only the results at 20 MHz have been discussed. The balance between refrac- 
tion and efficient growth is a function of frequency. At the highest frequencies, near 
27 MHz, even ^ s 80° is too small for convective growth because the very strong field 
gradients rapidly destroy the wave-particle resonance. At the other end of the spectrum, 
f 3E 1 0 MHz and below, the weak gradients in B permit substantial growth at all ^ ^ 50®. 
Below ^ ss 50®, our formula for y becomes invalid. 

This effect can be illustrated, if only approximately, by using Equations (14) and (15) 
to find the likely range of with which an amplified wave will leave the source region. 
The results at 27, 20, 16 and 10 MHz are shown in Figure 5. Values of 4^ larger than 
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about SS** are prohibited by refraction, while the lower limit at each frequency indicates 
where the convective amplification is small because changes in B have moved the particles 
out of resonance. 

The exact relationship between the observations and the results discussed above is 
not completely clear, for the observations sum over the entire source region. However 
the general trend of the theory is that at high frequencies (f ^ 20 MHz), radiation is pre- 
ferentially excited into a very thin sheet a 3®) with ^ a 85®; between 16 and 20 
MHz, ^ 10® and typically ^ a 80®. At still lower fr^uendes (f a 10 MHz) can 
range from 45-85®. In Paper II, we found that such a shift of 'if with frwjuency could 
ac(X)unt for the observed single peak in the occurren(» probability at 200® CML at 10 
MHz (Thieman and Smith 1978). 

These calculations could be improved by simultaneously solving Equations (10) and 
(15) along the ray path. The effects of refraction would then be rather predsely deter- 
mined, at least within the context of the field and density models employed. However, 
there are several complications with that approach, for as the ray is retracted, it moves 
away from the equatorial plane. Of course, some other ray not originally intersecting 
that plane will be refracted into it, but this means that to follow a ray toward a remote 
observer in the equatorial plane, one must actually move backwards along the ray path 
from the observer to the point where Km at that frequency is maximum. Such a calcu- 
lation is well beyond the scope of the present analysis, which we believe suffices in pro- 
viding reasonable estimates of refraction. 
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(b) Sourt^ Loations 

In this section we will use the results of Paper II along with estimates of K]^(X{|{) 
to find the predicted source location of DAM as seen by a distant observer in the Jovian 
equatorial plane. The numerical results discussed thus far have been evaluated at L = 6, 
the lo flux tube. We have also computed at L * 2, 3, 4 and 5. The range in 

over which Kj^i| is large does not vary with L, remaining between 80“ and 89,9“ (f > 10 
MHz). Furthermore, numerical estimates of both the amount of refraction and amplifica- 
tion of the ray again indicate that it will be refracted until ^ s 80“ , and that ^ ^ 80“ is 
precluded because gradients in B destroy the wave-particle resonance. Paradoxically, 
because these features of the analysis do not vary significantly with L, one can argue that 
the source region is, in fact, located near L = 6, First of all, by L ** 2 the growth rates 
have be<»me too small to account for the observed fluxes, but more importantly; we 
noted in Paper II that if the DAM source were inside L - 6, the geometrical model would 
better fit the locations of the Main and Early sources if ^ (or was no larger than 60“. 
However, as we have demonstrated above, at frequencies greater than 10 MHz, gradients 
in B destroy the resonance if { < 80°. Furthermore, we found that refraction is not 
likely to be strong enough to bend a ray by more than 10° at any L value. Thus, the 
results of the present analysis, taken together with the results of Paper II, suggest that 
the lo-independent sources are located near the lo flux tube. (Note that in neither Paper 
II, nor in the present analysis have L wilues greater than six been considered.) Exactly 
how lo might influence the excitation of the “lo-independent” DAM will not be consid- 
ered here. In the remainder of tliis subsection, we assume that the sources are located 
near L = 6 and use the O4 and JPL field models to determine the CML of the decameter 
sources. 
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In Figure 6 we show = 89*) at 20 MHz as a function of Xm uskg the 0^ 
magnetic field model. The results are shown for the northern hemisphere. For compari- 
son, 1 is also plotted. Note that the locations of large K|^ correspond to regions of small 
(. Thus we have a clear phyacal explanation for the relationship discusied in Paper II 
between t and source location, viz., in the northern hemisphere the loss-«)nc instability is 
strongest precisely where the axis of the emission a?ne has the greatest tilt toward the 
equatoral plane. 

From Figure 6a Kj^ is largest in the region 165® ^ Xm £ 210®. With ^ a 80®, 
radiation from this region will be beamed into the equatorial plane from one side of the 
(X>ne in the direction 

95® < A|j,(-/i) 170® 

From the other side of the cone radiation will be beamed toward 240® £ Ahj(+m) ^ 
315®. Over the part of the region where Am exceeds 270® the radiation is actually di- 
rected toward the cloud tops because of the large tilt of the cone axis out of the merid- 
ional plane. Therefore between Am(+*i) = 270®-315® radiation may be absorbed or 
reflected by the planet. We will make no allowance for this, and take the source region 
to be 

240® Am(+p) < 315® 

These ranges ^ould be compared with the observed locations of the non-Io B and A 
sources (Carr and Desch 1976 ~ this range refers to all frequencies between 1 1 and 28 
MHz, at 20 MHz the source sizes are smaller) 

95® A,i, ^ 195® Non-Io B 

195® < A„i 285® Non-Io A 
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In the southern hemisphere a similar calculation produces quite different results 
(Figure 6b). Now K|^ is relatively small where t is a minimum, and ampUflcation 
rates are foimd where the emi^on cone Ln tilted auay from the equat(»ial phne. The 
largest ^wth rates are found near X||j between 220^-260^. The resulting beaming pat- 
tmi, again a»uming ^ s 8(f , results in apparent source locations of 

290® % A|j|(+u) ^ 330® 
and 

155® % 195® 

while non-Io C is found between the (antisl meridian longitudes 

285® < A,„ % 5® 

From estimates of A* in these regions we found a tendency for A"*" to be significantly 
larger than A*. In Paper II, we showed that where t was small, the emission cone inter- 
sected the equatorial plane so that A|||(-u) corresponded to Source C. Now we find that 
the largest growth rates are found where the +u side of the cone corresponds to Source C. 
This is an encouraging development because it places the lo-independent source close to 
the location of the lo-controlled one. From the lo-phase versus CML relation^p, it is 
known that lo-C is emitted from the high longitude (i.e. +u) side of the emission cone. 
That is larger than A~ is consistent with the fact that no non-Io left polarized source 
is observed in the range 155-180® CML. 

Non-Io C is only observed at frequencies below about 23 MHz, At 23 MHz Ky has 
decreased substantially from its values at 20 MHz (Figure 6b). There is only one area 
(X||| ^ 235®) where Kj^ is not completely negligible. A^n A|||(+m) ^ 300® (Source C), 
and A'*' > A". Furthermore, in the region around Xjn =5 50®, where t is a minimum, 
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Kj^ ai 0. ThU tendr to confkm our associatiori of the emittii^ region of Source C with 
X|ll a 235^ in the iouthem hemiiphere. 

A timilar analy^ uting the JPL magnetic field model yiekts somewhat different re* 
suits. As noted in hiper II, this model diffo's substantially from in that the maximum 
southern hemisphere fields are found 100 ^ further west than the corresponding maximum 
in the O 4 model (d*. Figures 1 and 2 of Paper 10. This coupled with our basic a^mp- 
tion that the radiation is produced by trapped electrons rraults in rather different longi* 
tudM for the huge^ values of K]^. In the northern hemisphere (Figure 7a), Kj^ and t are 
no longer in phase. In fact, at the longitude of smallest t, is neat its minimum. 

In the northem hemisphere is largest between 220** % X|n iS 60** - a ran^ of 
some 20Cf. The r^uiting emisaon intta^n is beamed into the eqiatorial pkne from 

180® $ A|,i(-p) £ 300® 
and 

330® A,„(+^) ^ 80® 

where ^ » 80®. These results are not encouraging, for they do not approximate the 
known source locations of the Early and Main Sources. 

In the southern hemisphere Kj^ is large from Xju a 270® to X|n s 20®. Radiation 
beamed into an 80® cone will then be seen from the central meridian longitude ranges 

215® < Aiii(-p) < 310“ 
and 

350® ^ Ai|i(+/i) ^ 95® 

Although these do not appear to include the Source C region, we find evidence from e\^l- 
uating A^ that agnifiiant amplification does occur when A||j ^ 300-340®. Curiously, 
both A"*" and A" are large within that range. 
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Althou^ the O 4 Held mode! apparently provides a clearer explanation of the source 
locations, at least within the context of the assumption:: c-f this the<^, one should keep 
in miiKl that both mc^els are derive from similar data sets. In both cases the data sets 
are insufTicknt to predsely determine the Jovian magnetic field everywhere within the 
nuignetospl^re. In constructii^ the 0 ^ and JPL models, the limited nugnctic Held mea- 
surements were apparently analyzed using somewhat different analytical techniques. If 
one could asaibe “error bars” to the set of models derivable from the Pioneer 1 1 data, 
the O 4 and JPL models might then be representative ones - both equally valid to the 
predion of the data and limited spacecraft trajectory. Furthermore, our use of these 
models, especially in this paper, has tended to amplify the differences between them. 
Therefore one should probably take the view that it is encouraging that within the group 
of allowable ma:;netic field models, there appears to be at least one that is consistent with 
this theory of DAM. In addition, it is dear that detailed predictions of the source loca- 
tions of DAM sources are very sensitive to the local structure of the magnetic field . Only 
better measurements will resolve these differences, and allow closer comparisons Ixjtween 
theory and experiment. 

(c) Power Spectrum 

One aspect of the DAM observations not yet addressed is the way in which intensity 
varies with frequency. In Figure 8 the log of the peak intensity of DAM, F, is plotted 
apinst frequency. The curve, taken from Carr and Desch (1976), includes radiation from 
both lo-controllcd and lo-independent sources. However, in his thesis Desch (1976) 
found that when only the non-Io-conl rolled component is considered, the intensity spec" 
irum has the same shape, with reduced intensities. Only the shape of the curve concerns 
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As long as the dze of tl^ source regk>n is approximately (x>n$tant at all frequendes, 
the logarithm of the intendty should be roughly pro{K>rtional to the convective growth 
rate» Kj^. But Kj^ at a spedfic fr^uency is also a function of We have seen that at 
20 MHz K||(0 fe largest at { a 89.5® (Figure 3). As the frequency decreases so does the 
ai^le ^ at whi<'h is larged. This deorease is rather small until f falls below 1 0 MKz. 
From 10 to 2.5 MHz tl^ ^ of maximum decreases from 89® to 75®. Even smaller 
values are found below 2.5 MHz, but because the pwvth rate calculation assumes kjj/kj^ <5C 
I ({ ^ 45®) the resulting values for Kj^| are unreliable. The laigest K|^({) is also plotted 
in Figure 8 , the dot-da^ed line below 2.5 MHz indicating the frequency range in which 
{ is le^ than 75®. IThe calculation of Kj^| was done using the O 4 model along the north- 
ern hemisphere field line aossing the ^uatorial plane at Xm = 230®.] 

Because the relative scales between the two quantities were arbitrarily chosen, the 
agnificant features of the curves are that they both peak near the same frequency (f ss 
10 MHz for K|^, and f ^ 8 MHz for F), and thc^ both decrease to very small values at 1 
and 27 MHz. 

The shape of the Kj^ curve is primr. .y dictat. 4 by the assumption that the instability 
is exdted by free energy in the particle’s loss-cone distrib 'ion. At the lowest frequen- 
cies, the lotal loss-cone angle is small aj 15°) so there is little velocity-space aniso- 
tropy to drive the instability. Conversely, at the highest froqucncies, is large ^ 

75° at 27 MHz) but nearly all electrons have mirrored, leaving too few particles to excite 
a strong instability. Somewhere in between these two competing effects convolve to pro- 
duce the maximum instability. The shape of the observed peak tlux curve is evidence 
that the particles’ loss-cone angle plays ; n important role in exciting DAM. 



(d) Emitting Region and Intensity 


Jovian noise storms have been observed using very long l^seline interferometer 
(VLB!) radio antenna systems in an attempt to measure the extent of the emitting repon. 
Obsc»‘vations have been confined to looking at the lo-B source at 34 and 18 MHz. At 
34 MHz, the observed noise storms were S-bursts (Dulk 1970), while Lynch et al. (1972), 
observing at 18 MHz, analyzed an L-burst. In both cases, the sources, assuming they 
were incoherent, were unresolved, indicating linear dimensions less than 400 km. To our 
knowledge, no interferometry has been done on non-lo sources. For the purposes of this 
discu^ion, tl^n, we assume that the non-lo sources al«) have linear dimen^ns of several 
hundred kilometers. Within such a distance the radiation must be amplified from noise 
to a peak flux of nearly 10^ Jy* (at 20 MHz). That noise level can be estimated by as- 
suming that electrons participating in the instability amplify their own incoherent cyclo- 
tron radiation. At earth, the noise level would be (Jackson 1962) 
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where (Ax)“ is the area of the miitting region (taken to be 1.6 x 10^ km"), Af is the 
receiver bandwidth and 
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With Ax = (400 km), F^^ 1 Jy, which implies that A' must be about 16. Such ampli- 

fications are rather easily adiicved at the longitudes characteristic of the peak values of 
Km- Furthermorwrthese amplifications typically result from traversing a path length of 
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a few hundred kilometo^. Thus the growth rates we have computed appear fully caj^ble 
of amplifying incoherent cyclotron taxation to the peak observed intensities within di- 
men^ons le^ than, or of order, the upper limits of tlte source size. The time necessary 
to reach the^ intensities is less than one %oond. 

That the excited waves convect out of resonan<^ after ampUfyii^ by the approximate 
amount required to account for the ob^rv^ intensities su^ests that nonlinear processes 
may not always be important in limiting the growth of the waves. Also, becau^ we lack 
sufficiently detailed information about the plasma environment in the emitting region, it 
is difficult to suggest mechanisms whereby nonlinear processes, be they wave-i»rticle or 
vrave-wave, limit the growth of the instability. 

There is, however, one effect that will tend to inhibit amplifications much in excess 
of the e*^ estimated above. Because the growth rates are so sharply peaked in frequency, 
electrons will become phase wherent or magnetically trapped by the wave if the ampli- 
tude becomes too large. Equations (16) and (17) can be used to estimate the noise level 
of the magnetic field fluctuations. At 20 MHz we find that 8B as 10"* ^ G, and the mag- 
netic trapping time, t is of the order 
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With 7 ~ lO'*'^ s"* there is time for 10 or so exponentiations before t exceeds 7"* , tend- 
ing to inhibit further wave growth. 


IV. SUMMARY 

The theory presented here has evolved from one basic assumption; the inner Jovian 
magnetosphere (L 6) is populated by keV electrons stably trapped in Jupiter’s magnetic 
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ftekt. From this, we have shown that the electron distribution is unstable to the elec^o- 
magnetic lo^-cone instability, producing radiation in the R-X mode near the local elec- 
tron gyrofrequency. Significant amplification is expected from 1 MHz to 28 MHz, with 
a peak near 10 MHz, which corresponds very closely to the observed frequency range and 
spectral shape of the lo-independent decameter spectrum. The competition between 
refraction and convective growth of the radiation in regions of strong radial gradients in B 
limits the radiation pattern at frequencies above 10 MHz to a wide, but thin, conical sheet 
with ^ = 80®-85®. At altitudes where tlw gymfr^uency is le^ than or ^ual to 10 MHz, 
the gradients in B are relatively small so that significant wave growth is possible over a 
much wider range of cone angles: 50® 'k ^ 80®. 

Thieman and Smith (1978), showed that by 10 MHz sources A and B appear to have 
been replaced by a single source near Am 2 200®. We argued in Paper II th.^t the effect 
would result from a decrease in 'P from 80° above 10 MHz to 'P a 40® at 10 MHz. Effi- 
cient wave excitation at frequencies below 10 MHz with »P as small as 50® or less is cer- 
tainly suggested by our results. 

Because of the very non-dipolar nature of the Jovian magnetic field, large variations 
of K}^ with Xm were found which enabled us to predict the location of the DAM sources. 
Those results are summarized in Table 1. Radiation excited where Kjyj is largest is beam- 
ed into the equatorial plane from either side of the wide emission cone ('P = 80°), From 
the northern hemisphere (O 4 model), the apparent source locations agree quite well with 
the known locations of non~Io A and B to within 30°. Using the JPL model results in a 
poorer fit. From the southern hemisphere -both field models suggest a DAM source con- 
sistent with the location of non~Io C. In addition the model indicates the possible 
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exij^ence of a weak left polarized sour^ near 155°- 1 80° CML. Such a source, if it exists, 
should be observable by the Voyager spacecxaft. Furthermore we found tlat the theory 
requires that the lo-independent DAM be excited on flux tubes with L ^ 6. 

In Paper II we noted that the geometry of the northern hemisphere sources could be 
understood by associating the source locations with small values of i. This postulate has 
a natural explanation when the field model is used, for then i is smallest where Kj^| is 
largest (Figure 6a), That the postulate failed to account for the location of source C can 
be understood because in the southern hemisphere small values of C are no longer corre- 
lated with large Kj^|, and the source location for non-Io is more naturally associated with 
the large values of Kj^| (Figure 6b). 

Because the time scales for exponential growth are of the order 10'^ s, S-bursts 
cannot be excited by this mechanism, but S-bursts are a unique feature of the lo-con- 
trolled sources; so this is not a limitation of the theory. The growth rates we have found 
are also consistent with constraints placed on the source size by VLBI observations, as- 
suming those results can be extrapolated to the lo-independent sources; i.e., the observed 
peak fluxes are achieved after convective growth over distances of 400 kn or less. At 
20 MHz, larger source sizes are unlikely because the gradients in B wilt move a wave out 
of resonance by the time it has travelled 400 km. 

(a) Some Unsolved Problems 

The theory presented here appears capable of explaining many features of decameter 
phenomenology. However, there are aspects of the observations that we have not addressed. 
Chief among them is the well known asymmetry in the occurrence probabilities of the 
non-Io A and B sources. Even in the Desch et al. (1975) study at 26.3 MHz, where 
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non-Io B had an occurrence probability twi<^ that of non-lo A, the intensities of the two 
sources differed greatly (B was less intense than A). The explanation of this effect may re- 
quire a (x>mbination of better magnetic field models and more precise ray-tracing calculations 
than used here. Our estimates of A* (Eq. 1 5) did suggest differences between A"^ and A", 
but none sufficiently striking as to provide a clear explanation for the observed differences 
between non-Io A and B. 

Another problem not discussed is the question of the source of the trapped electrons. 
The instability we (x>nsidered will stabilize and turn itself off at any particular location 
after several seconds. Additional bursts from that area would require a fresh supply of 
electrons. If the L = 6 flux tube is the source of the radiation, then lo could perhaps 
indirectly provide the necessary particles. Alternatively, the suggestion has been made 
that the occurrence of lo-independent emission is correlated with fluctuations in solar 
wind parameters (v., e.g. Terasawa, Maezawa, and Machida 1978), suggesting the solar 
wind as the possible source of the electrons. This brings us to the question of the rela- 
tionship between this mechanism for exciting the lo-independent DAM and theories of 
the lo-controlled component. It is a question we hope to address in future research. 
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Table 1 


Comparison of Observed and Predicted Source Locations Uaiiig Both the O 4 
and JPL Multipole Magnetic Field Models 


Observed Source 
Locations (after 
Carr and Desch 
1976) 

Polarization 
R = Right, L = Left 

f»20 MHz,4f-80°,L«6 

O 4 

JPL 

Observed 

Predicted 

Source Locations 

^III of 

Emitting Region 

Source Locations 

hu of 

Emitting Region 

Non-Io A 







195°-285° 

R 

R 

Ain(+p) 


Aiii(+|i) 


(11-28 MHz) 



240°-315° 


330°-80° 






165°-210° 


22<r-60“ 

Non-Io B 







95°-195° 

R 

R 

Aui(-m) 


Aiii(-*t) 


(11-26 MHz) 



95°-170° 


180°-300° 


Non-Io C 







285°-5° 

L 

L 

Aiii(+m) 


Aiii(+M) 


(11-22 MHz) 



290°-330° 


350°-95°t 






200°-260° 


270°-20° 



L 

Aiii(-m) 


Aiii(-tM) 





155°-195°* 


215°-310°t 



*Estimates of A'* indicate weak amplification. 

tEstimates of and A~ indicate large amplification in the 300-340° longitude range. 
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Figure 1 . Relationship between Vg, (o/k and Note that 
when the Appleton-Hartree dispersion relation is used, 5 is 
negative for the extraordin’^ry mode in Band III. 
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Fig»e 2. C«>w,h «.c f«j,u.ncy ta M 

^ ■ .iL o!, mMel MThTpoin. a, = 56" oomp«r«J to «, - 2.5* a. ttc n»g»eUc equator. 





Figure 3 . Variation of K„ with J for several values of the power-law index, m. N, { and 
® Ul as well as r, 0, and Xm are the same as in Figure 2. 
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Figure 5. Estimated variation of with frequency. The limit for large angles is determined by refraction. The cutoff 
toward small angles indicates when convective amplification is small because the wave has moved out of resonance 

with the electrons. 
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O4 NORTHERN HEMISPHERE 




Figure 6 . A plot of Kj^| and t versus Xm in the northern (a) and southern (b) 
hemispheres using the 64 magnetic field model. Note that K]^ is large in the 
northern hemisphere where i is small. In the southern hemisphere Kj^ is large 
where t is large. Note that in the southern hemisphere t is plotted for the 
20 MHz source height, while is given for both 20 and 23 MHz, 
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JPL NORTHERN HEMISPHERE 


40 ' 



JPL SOUTHERN HEMISPHERE 



7. Similar to Figure 6, but using the JPL magnetic field model. Note that 
i and K[yi are not well correlated in either hemisphere. 
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Figure 8. Peak radio flux and Kj^ plotted against frequency. The radio data is taken from 
Carr and Desch (1976). Kj^ is the largest value of at each frequency. The shape of 
the Km curve is a consequence of using a loss-cone distribution to excite DAM. 
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TABLES 


1 Comparison of observed and predicted source locations using both the O4 and JPL 
inultipole magnetic field models. 


ILLUSTRATIONS 

1 Relationship between v^, <j/k and B Note that when the Appleton-Hartree disper- 
sion relation is used, 6 is negative for the extraordinary mode in Band III. 

2 Growth rate versus frequency in Band III. In this calculation m = 2, N = 100, f = 
0.1 and u^ = 5.9 x 10* cm 8"^. B was evaluated at r = 1.06 R^^, 0 = 54“ north 
latitude, and Xjn = 201.7“ w iih the O4 model. At that point = 56” compared to 
ot£ = 2.5“ at the magnetic quator. 

3 Variation of with { for several values of the power-law index, m. N, f and ul 
as well as r, d, and X{]| are the same as in Figure 2. 

4 Variation of TM/lfle ^ ^ several values of m. Other parameters are the same 

as in Figure 2. 

5 Estimated variation of 4^ with frequency. The limit for large angles is determined 
by refraction. The cutoff toward small angles indicates when convective amplifica- 
tion is small because the wave has moved out of resonance with the electrons. 

6 A plot of K^i and i versus X|j| in the northern (a) and southern (b) hemispheres 
using the O4 magnetic field model. Note that Kj^ is large in the northern hemi- 
sphere where i is small. In the southern hemisphere Kj^ is large where i is large. 
Note that in the southern hemisphere i is plotted for the 20 MHz source height, 
while Km is given for both 20 and 23 MHz. 

7 Similar to Figure 6, but using the JPL magnetic field model. Note that i and K^ 
are not well correlated in cither hemisphere. 

8 Peak radio flux and plotted against frequency. The radio data is taken from 
Carr and Desch (1976). K^ is the largest value of K^ft) at each frequency. The 
shape of the curve is a consequence of using a loss-cone distribution to excite 
DAM. 
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